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ABSTRACT 

We describe a method to measure the A^bh-c* relation in the non-local universe using dust- 
obscured QSOs. We present results from a pilot sample of nine 2MASS red QSOs with redshifts 
0.14 < z < 0.37. We find that there is an offset (0.8 dex, on average) between the position of 
our objects and the local relation for AGN, in the sense that the majority of red QSO hosts 
have lower velocity dispersions and/or more massive BHs than local galaxies. These results are 
in agreement with recent studies of AGN at similar and higher redshifts. This could indicate an 
unusually rapid growth in the host galaxies since z ~ 0.2, if these objects were to land in the local 
relation at present time. However, the z > 0.1 AGN (including our sample and those of previous 
studies) have significantly higher A^bh than those of local AGN, so a direct comparison is not 
straightforward. Further, using several samples of local and higher-z AGN, we find a striking 
trend of an increasing offset with respect to the local A^bh-c* relation as a function of AGN 
luminosity with virtually all objects with log(L5ioo/erg s~^) > 43.6 falling above the relation. 
Given the relatively small number of AGN at z > 0.1 for which there are direct measurements 
of stellar velocity dispersions, it is impossible at present to determine whether there truly is 
evolution in A^bh-c* with redshift. Larger, carefully selected samples of AGN are necessary to 
disentangle the dependence of A^bh-c* on mass, luminosity, accretion rates, and redshift. 

Subject headings: galaxies: active — galaxies: evolution — quasars: general — dust: extinction 



1. Introduction 

During the last decade, the astronomy com- 
munity has come to realize that black hole (BH) 
activity plays an important role in galaxy forma- 
tion, e.g., by regulatin g star formation t h rough 
winds and outflow s ( Scannapieco et all l2005t 
Cattaneo et al.|[2005l) . To understand how galax- 
ies form, we therefore need to understand how 
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supermassive BHs and galaxies co-evolve. 

The most important observational tools for 
studying the co-evolution of BHs and galaxies 
are the empirical scaling relations between BH 
mass and global ga laxy properties, such_ a.s the 



relation (Ferrarese and Merritt 2000l: 



Gebhardt et al.1 l2000al) . a particularly tight re- 
lation between BH mass, AIbh, and the cen- 
tral, stellar velocity dispersion, a.^,. The A^bh- 
cr* relation suggests a causal connection be- 
tween the BH and the b ulge (although see e.g., 
Jahnke and Macci6ll2011 . for a counter argument). 



and is believed to hold clues to understanding the 
galaxy formation process. By studying how the 
BH mass relates to the mass of the bulge (as mea- 
sured by cr*), the galaxy evolution process can be 
examined across different redshift ranges and for 
different galaxy types. 

Active galactic nuclei (AGN) have become in- 
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strumental in this work because they contain ac- 
tively accreting black holes. Not only are they 
direct probes of the BH-galaxy co-evolution, but 
galaxies with AGN are the only galaxies in which 
we can measure black hole masses in the non- 
local universe. The method is referred to as the 
virial method, where the width of broad AGN 
emission lines s uch as Ha or H/3 is used to mea- 
sure A^BH (e.g.. IVestergaard and Peters on 2006). 
Whereas measuring virial BH masses in AGN 
through broad emission lines is relatively straight- 
forward, measuring the host galaxy velocity dis- 
persion is hampered by the presence of the very 
bright nuclei which can completely overpower the 
light from the underlying galaxy. 

Attempts have been made to use the width 



of [Q UI] line s as a proxy for the bulg e potential 
dNelson 2000; IShields et al.. .2003: Salviander et al 



2nn7f i. 



but these studies suffer from a large 
scatter due to non-gravitational influences on 
the [OIII] line width (such as outflows, ra- 
dio luminosity, accretion rate s, etc; see, e.g., 
iNetzer and Trakhtenbro^ l2007t ). Studies of low- 
luminosity AGN at z < 0.1 have shown that 
they follow the A^bh-c* relation of loca l , more 

2000bt 



massive spheroids 



...^...v^ ^^..^.^.^^ (e.g., — ■ 

Nelson et al. 20041: iGreene a nd Ho 2006b HWoo et al 



Gebhardt et alj 



20101 : iBennert et al.ll2011al ). and there is evidence 
that higher redshift AGN show evolution with 
regard to the local TWeh-c* relation. 



Woo et all (|2006l l2008l) and iTreu et all (|2007l ) 



studied two samples of Seyfert 1 galaxies at z = 
0.36 and z = 0.57 and found an offset from the 
local relation suggesting that, for a fixed black 
hole mass, higher redshift spheroids have smaller 
velocity dispersions than local ones. This indi- 
cates that, in these objects, the gala ctic bulge is 



still f orming around the black hole. iHiner et al 
(j2012l) find a similar offset with respect to the lo- 
cal relation in a sample of six post-starburst AGN 
at z ^ 0.3. These results have serious implica- 
tions and may provide important observational 
constraints for galactic evolution models which at- 
tempt to explain the co-e volution of black holes 
and bulges. The studies by lWoo et al. . Treu et al 



and IHiner et al.l suggest a very recent growth, 
within the last 5—6 billion years, of intermediate- 
mass bulges for a given black hole mass, in con- 
trast to predictions ma de by detailed hydrody - 
namic simulations (e.g.. [Robertson et al.l l2006bl ) . 



On the other hand, IShen et al.l (|2008l ) study a 
large sample of AGN observed with the Sloan Dig- 
ital Sky Survey (SDSS) at 2 < 0.452 and find no 
evidence of evolution in the relation. 

The study of the A1bh-c* relation has been 
severely limited in QSOs due to the observational 
challenges such a project poses. However, study- 
ing A^bh-c* in QSOs would be desirable since 
these are potentially the objects where the most 
rapid growth (in both BH and bulge) occurs. 
Studies of A^bh-o"* at z > need to include QSOs 
in order to probe the higher end of the relation 
and, potentially, episodes when either the bulge 
or the BH is catching up in its growth. 

Measuring velocity dispersions from absorp- 
tion lines in QSO host galaxies is observation- 
ally challenging. The relatively few studies of 
host galaxy stellar a bsorption line spectroscopy in 



the literature (e.g ., Canalizo and StocktonI 200lL 



'N olan et ahlbOOirijahnke et al.l 120071: IWold eTal 



|2010[ ) are not able to recover host spectra in the 
central regions of the galaxy. Using these spec- 
tra to measure velocity dispersions could lead to 
ambiguous results since it is difficult to determine 
whether the dispersions are being measured in the 
bulge of the galaxy or in the outskirts or some ex- 
tended tidal feature where cr* may not be represen- 
tative of the mass of the bulge or that of the BH. 
Indeed, few attempts have been made to measure 
a^: in the host galaxies of luminous QSOs (e.g.. 
Wolf and SheinislbOOsHRothberg et al.|[2012l ). but 



these studies have been in relation to the host 
structural properties rather than A^bh-c* . There- 
fore, it would be ideal to find a class of objects 
with a natural "coronagraph" that would occult 
the QSO nucleus when measuring cr*, but one that 
we could remove in order to measure A^bh from 
the width of the broad lines. As described below, 
we have found that dust-reddened QSOs fit the 
biU. 

Red QSOs have some of the same characteris- 
tics as their blue counterparts, such as (at least 
some of the) strong broad emission lines and high 
bolometric luminosities, but with much redder 
continua. At present, there is not a clear defini- 
tion for red QSOs, so that the different objects 
that are cataloged as red QSOs do not form a 
homogeneous class. The majority of red QSOs 
discovered to date, however, appear to be red- 
dened by dust, and thus they are considered the 
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lation fe.K., Cut 


ri et al.ll2002l; 'Marble et al.| 


20031 


Hall et al. 20021 


Glikman et al. 2004: White et al. 


2003). Even then, there are cases where the dust is 



intervening rather than intrinsic to the QSO (e.g., 



Gregg et al.ll2002l ). Here we focus on those objects 
where the dust is presumably near the nucleus. 

The heavy extinction of the QSO nuclei at opti- 
cal wavelengths makes dust reddened QSOs excel- 
lent candidates to study velocity dispersions and 
stellar populations in their host galaxies. The 
host galaxy spectrum suffers from much lower (if 
any) extinction than the nucleus, and stellar ab- 
sorption features such as those around the Mgib 
and G-band regions are clearly visible in the blue 
region, being almost free of contamination from 
QSO light. At longer wavelengths, however, the 
contribution from the reddened QSO continuum 
increases and broad Ha is clearly visible. Thus, 
cr* can be measured from features in the underly- 
ing host galaxy and A^bh can be estimated from 
the width of Ha using the same spectrum. 

In this paper, we present a pilot study to test 
the feasibility of using these objects to measure 
A^BH-c*. In § [5] we describe the sample, imag- 
ing and spectroscopic observations, as well as the 
data reduction. In § [3] we describe our method 
to fit stellar velocity dispersions and to measure 
A^BH- In §|4]we present our results and consider 
potential biases in our sample, while in § [5] we 
discuss our results in the context of other studies 
of A^BH-c* that use AGN. Throughout the pa- 
per, we have assumed a cosmological model with 
■ ffn = 71 kms~/Mpc ~\ VL^ = 0.27 and VLa = 0.73 



(jSpergel et al.ll2003[ ). except where noted. We use 
the Vega magnitude system unless otherwise spec- 
ified. 

2. Observations and Data Reduction 
2.1. The Sample 

Our sample is drawn fro m the sample of 29 
red 2MASS QSOs studied bv lMarble etaD (|2003), 
w hich in turn i s draw n from the sample compiled 
bv lCutri et al](l2002h. We focused on the subsam- 



ple from iMarble et al.l (|2003[ ) because these ob- 
jects have Hubble Space Telescope (HST) imag- 
ing observations that are useful to measure the 
nuclear luminosities and to study the properties 
of the host galaxies. The objects in this sam- 



ple have Mk < -25, J ~ Kg > 2.0 and detec- 
tions in each of the three 2MASS bands, JHKg- 
They have a median redshift of z = 0.213 and 
have beei i spectroscopically confirmed to contain 
an AGN (|Smith et al.ll2002l ). 

A fraction of these objects hav e spec tropo- 
larimetry published bv ISmith et al. (2003). We 
avoided the objects that show a blue continuum, 
with broad emission features dominating the op- 
tical spectrum. Our final sample consists of nine 
of the remaining targets chosen at random within 
the observational constraints. 

In Table [1] we list the final sample. The red- 
shifts listed were measured from stellar absorp- 
tion lines in our Keck spectra (§ 12. 2p . J — Kg 
color s and absolute Ks magnitudes are quoted 
from iMarble et al.l ( 2003 ) , and the latter assume 
Hq = 75 kms"^Mpc"^ and go = 0.5. 

2.2. Spectroscopy 

We obtained deep, medium resolution spectro- 
scopic observations wit h the Echellette S pectro- 
graph and Imager fESl dSheinis et al]|2002[ ) on the 
Keck II telescope during three separate nights in 
October 2003, July 2004 and September 2004. We 
observed in the echellette mode, which provides a 
wavelength coverage from 3900 to 11000 A in 10 
orders (6 through 15). We used a 1" slit, which 
projected to ~7 pixels (ranging from ^^8 pixels 
at the short-wavelength end to '^6 pixels at the 
long-wavelength end) on the MIT-Lincoln Labs 
2048x4096 CCD detector. The spectral resolution 
is 11.4 km pixel"^. We placed the slit through 
the center of the host galaxies to ensure that the 
velocity dispersions that we measure correspond 
to those of the host galaxy bulges. 

All targets were observed under clear weather 
conditions and subarcsecond seeing (~ 0'.'6 in V), 
except for 230442.4-^270616 and 022150.6-^132741, 
which were observed through clouds, and 232745.6 
-1-162434, which was observed in high wind condi- 
tions. Exposure times for each target ranged from 
1800 to 7200 s; specific times are hsted in Table[T] 

We also observed a suite of template stars with 
stellar types ranging from FO to M6 (listed on Ta- 
ble [2]) and several spectrophotometric standards. 

The spectra were reduced using a combination 
of IDL and IRAF scripts written to meet the spe- 
cific needs of our program. After bias subtraction 
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Table 1 
The Red QSO Sample 



Object 


QSO 


Zhost 






ESI Exposure 


Extraction aperture 


ID 


(2MASSi J) 








Time (sec) 


(arcsec) 


(kpc) 


1 


005055.7+293328 


0.1356 


-25.6 


2.1 


1800 


0.21 


0.50 


2 


015721.0+171248 


0.2139 


-27.0 


2.7 


7200 


0.21 


0.75 


3 


022150.6+132741 


0.1398 


-25.7 


2.4 


5400 


2.86 


6.96 


4 


034857.6+125547 


0.2112 


-26.7 


3.3 


3600 


0.79 


2.70 


5 


163736.5+254302 


0.2772 


-26.5 


2.3 


5400 


0.21 


0.88 


6 


165939.7+183436 


0.1709 


-26.5 


2.2 


5400 


0.29 


0.84 


7 


225902.5+124646 


0.1989 


-25.6 


1.9 


5400 


0.42 


1.36 


8 


230442.4+270616 


0.2370 


-25.4 


2.1 


7200 


0.97 


3.61 


9 


232745.6+162434 


0.3664 


-27.0 


2.4 


5400 


1.29 


6.51 



Table 2 

Stellar templates used for the velocity dispersion fitting 



Star 


Luminosity 




class 


HD 218140 


FO 


SAO 90936 


F5 


HD 11851 


F8 


BD+264556 


GO 


HD 10995 


G2IV 


SAO 91026 


G5 V 


SAO 090989 


KO 


SAO 91028 


K2 


HD 11326 


K2III 


HD 218113 


K5ni 


SAO 90990 


K5 V 


0120131 


MO 


HD 11729 


M6 


SAO 92712 


KO 


SAO 92718 


GO 
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and flat fielding, each order was rectified indepen- 
dently in both spatial and wavelength directions 
using a modified version of the WMKONSPEC 
packag^ll- A star was observed at multiple po- 
sitions across the slit in order to trace and rec- 
tify the spatial direction of the spectra. A wave- 
length solution, obtained from observations of a 
Cu-Ar lamp, was used to rectify the wavelength 
direction. Sky lines were removed by fitting the 
background in the two-dimensional spectra. The 
spectra were c alibrated using spectrophotometric 
standards from lMassev et al.l ([1988) observed with 
the slit at the parallactic angle. Each target had 
two or three individual exposures; we averaged the 
spatially corrected spectra using the IRAF task 
scombine. 

We extracted spectra from the central few kpc 
of each spectrum, adjusting the size of the aper- 
ture to match approximately the effective radius, 
Tcff, of each target (see § 13.21) . as fisted in Ta- 
ble [H This corresponded to a slightly different 
aperture size for each echellette order, since the 
spatial scale varies for different orders. We then 
combined the spectra from all orders into a single 
spectrum for each target. We corrected the spec- 
tra for_G^la£feextmction, using the values given 
bv ISchlegel et al. ( 19981 ). Finally, we transformed 
the spectra to rest frame using redshifts measured 
from stellar absorption lines. The final spectra for 
the nine targets are shown in Figs. Hal and llbl 

2.3. Imaging 

The imaging data were obtained from the HST 
data archive (SNAP-90 57; PI: D . Hincs) and were 
originallv published bv lMarble et al.i (|2003h . The 
targets were observed with the PC chip of the 
WFPC2 in the F814W filter. Each target had two 
400 s exposures. We combined the images with 
multidrizzle, using standard procedures. 

3. Analysis 

3.1. Velocity Dispersion 

The velocity dispersion in the host galaxies 
was measured by fitting broadened stellar tem- 
plates to the galaxy spectra. A dire ct fitting al- 
gorithm based on that described by iBarth et~al 




Fig. la. — Keck ESI rest-frame spectra of our sam- 
ple of nine 2MASS QSOs. 



()2002l) was used. The algorithm has been exten 
sively te sted and utilize d for A GN h ost galaxies 
by, e. g., iGreene and Ho (|2nn6a^ and IWold et al 
(|2007l ). The model spectrum M[x) which is fit- 



-"^http: / /www2.keck.hawaii.cdu/inst/nirspec- 
old / wmkonspcc /index. html 



ted to the host galaxy is formed by convolving 
a stellar template T{x) with a line-of-sight veloc- 
ity broadening function, assumed to be a Gaussian 
with dispersion , adding a power- law continuum 
C{x) and multiplying by a polynomial P{x): 

M{x) = {[T{x) ® G{x)] + C{x)} X P{x), (1) 

where T{x) ® G{x) is the stellar template con- 
volved with the Gaussian, G{x). The x coordinate 
is defined as a; = In A so that velocity shifts are lin- 
ear in In A. The continuum is C{x) = cq -I- cix 
(note that this is a power-law as a function of 
wavelength), and P{x) was chosen to be a 3rd 
degree polynomial P{x) = C2Pq(x) -I- c^pi{x) + 
C4P2 (x) + C5P3 (x) , where pn (x) are Legendre poly- 
nomials. In order to find the best-fitting model, 
we varied the velocity dispersion and the six coef- 
ficients cn to c.g using a downhill simplex method 
( Press and SpergelllQSSi) until a global minimum 
of the function was found. 

Our object spectra are combinations of QSO 
and galaxy features, and the featureless contin- 
uum C{x) describes the underlying AGN contin- 
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Fig. lb. — Keck ESI rest-frame spectra of our sam- 
ple of nine 2MASS QSOs. 



uum. The polynomial P{x) accounts for reddening 
between the template and the object, so the factor 
C{x) X P{x) which is added to the convolved stellar 
template describes the underlying reddened QSO 
spectrum. By comparing C{x) x P{x) with an un- 
obscured QSO spectrum, we were able to obtain 
an estimate of the amount of reddening for each 
QSO. This is described in §[3X2l 

Stellar templates T{x) were formed from indi- 
vidual and various combinations of the 15 different 
stellar spectra observed with the same setup as the 
red QSOs (see Table [H). Average F, G, K and M 
spectra were made by normalizing to the flux in 
the 7000-7100 A region and taking the average of 
all spectra within each of the four classes. Since 
the main source of uncertainty in measuring ve- 
locity dispers ions typically come s from template 



misma tches (jBarth et al.l 120021 : iGreene and Ho 
l2006ah . we ran a series of tests to determine the 
best template to fit for each region and thus min- 
imize systematic uncertainties, as described in 
detail below. 

Based on previous experience and descriptions 



of similar analyses in the literature ([Barth et aL. 
20021 : iGreene and Holl2006al : IWold et al.ll2007h we 



chose to estimate velocity dispersions from fits to 
two different spectral regions: (1) a blue region 
at 3900-4600 A including CaK but excluding the 
CaH line (the region 3950-3985 A around CaH 
was excluded because in some cases AGN con- 



tamination from He in emission complicated the 
fit), and (2) a red region at 5100-5550 A with 
the Mg lb absorption line complex AA5167, 5173, 
5184 masked out (5155-5220 A). It is well-known 
that the correlation between [Mg /Fe] and velocity 
dispers ion for elliptical galax ies (|j0rgensen et al "" 



[19991 : iKuntschner et al.ll200l[ ) can cause problems 
with simultaneous fitting of th e Mg lb absorption 



and the region r edward of it (jBarth et al. 1 120021: 



Wold et al.l l2007f ) . The red region may also be 
affected by the AGN Fe II pseudo-continuum at 
5050—5520 A, affecting the fits to this region. This 
is seen in a couple of our objects, as well as emis- 
sion from [N I] AA 5198, 5200. We therefore chose 
to mask out the region 5155-5220 A from our fits. 
Apart from this, the chosen red and blue regions 
are the two most reliable regions in the optical 
with strong enough stellar features to be used for 
velocity dispersion measurements. 

The region redward of the Mg lb contains sev- 
eral iron absorption line features that are well fit- 
ted with our composite template. We also obtain 
good fits to the blue region, in particular to the 
CaK line at 3933 A, H5-absorption at 4100 A, and 
the G-band at 4300 A. Generally, we do not see 
much AGN contamination in this region, except 
in one case (165939-f 183436), which consequently 
has larger uncertainties in the velocity dispersion. 
Some narrow AGN emission lines, such as H7, H5, 
and [O III] A4363, had to be masked from some of 
the fits, as well as corrupt spectral regions in a cou- 
ple of cases. In cases where we could not achieve a 
good fit to the 4OOOA break, but there were enough 
features longward of the break to constrain the fit, 
we restricted our fit to a region from 4000—4600 

A. 

The galaxy spectra were shifted to the rest 
frame of the stellar templates, and the fitting algo- 
rithm initiated with a reasonable set of parameters 
(we fit seven free parameters including the veloc- 
ity dispersion) . The parameters were varied freely 
using the downhill simplex routine until a best fit 
was found. To improve on the best fit and to en- 
sure that a global minimum was found, the fitting 
routine was started again, with initial model pa- 
rameters equal to the first solution but perturbed 
by 10—20 % in random directions for each param- 
eter. The confidence interval on was thereafter 
found by varying the velocity dispersion in steps of 
10 km s^^ while letting the other six parameters 
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float freely. 

In order to minimize systematic uncertainties 
due to template mismatches, we used single stars 
as well as different combinations as stellar tem- 
plates. For the red region, we find that, in gen- 
eral, we obtain the lowest reduced when using 
late stellar types as templates. Five of the ob- 
jects were best fit in this region with a K 5 tem- 
plate. This is not surprising since this region is 
generally dominated by features from later stel- 
lar types. This also consistent with the results 
found by Hiner et al. (in preparation), who mea- 
sure stellar velocity dispersions in local galaxies by 
using templates constructed from a large suite of 
stellar types, accurately mimic king stellar popula - 
tions of different ages (see also lHiner et al.ll2012l) . 
They find a one-to-one correlation in the veloc- 
ity dispersions that they measure by using a sin- 
gle K star template versus more complex stel- 
lar templates in the Mg lb region (equivalent to 
our "red region"). Only three of the objects in 
our sample (005055. 7-H293328, 034857.6+125547, 
225902.5-H124646) are best fit by including a con- 
tribution from earlier stellar types. Although 
these fits produced the smallest reduced x^: the 
resulting cr* is within a few tens of km s~^ of 
that measured with a K star template. In fact, 
the standard deviation in the values of cr* mea- 
sured from different templates is typically smaller 
than ^20 km s~^. Since we do not do a detailed 
analysis of the stellar populations in this work, we 
do not favor any particular composite template 
among those that yield the lowest reduced val- 
ues. Therefore, for consistency, we chose to use the 
same composite template for these three objects as 
the one we used for the blue region, as described 
below. The remaining object, 015721.0-1-171248, 
shows significant QSO contamination in the red 
region and we were not able to obtain a reliable 
fit in this region. 

The spectra of our objects in the blue region 
show features from both early- and late- type stars. 
Therefore, fits using single stellar types as tem- 
plates consistently resulted in significantly larger 
reduced ^ values than those obtained from fits 
using composite stellar templates. This is also 
consistent with the results by Hiner et al., who 
find that using single stellar types as templates 
for this spectral region result in systematic uncer- 
tainties of several tens of km s^^. Thus, we tested 



a few templates consisting of combinations of stel- 
lar types with ratios by flux that would mimic 
intermediate-age populations (i.e., including fea- 
tures from older and younger populations) and 
chose the one that yield the lowest values of re- 
duced x^- Although, in principle, we could do a 
X^ minimization to determine the best combina- 
tion of stellar types for each template, we chose 
use the same template for every host galaxy in or- 
der to reduce the number of free parameters (and 
degeneracies) in the fits. The flux contribution 
at 7000 A of each stellar type for this composite 
template is 40% K, 40% G, 15% F and 5% M-star 
(Fig. El. 




Fig. 2. — Composite stellar template, consisting 
of 40% K, 40% G, 15% F and 5% M-type stars 
(by flux). The stellar template shown here has 
been convolved with a cr = 165 km s^^, which is a 
typical velocity dispersion of the host galaxies in 
the sample. 



The galaxy spectra in the red and the blue 
regions with the best-fit models overplotted are 
shown in Figs. l3al and I3b| and the velocity disper- 
sions obtained from these fits are listed in Tabled 
The reduced chi-squared, Xred' fo'" each fit is also 
listed. 

The confidence intervals on cr* give the un- 
certainty related to the fitting procedure, with 
smaller confidence intervals for higher signal-to- 
noise spectra and for spectra with minimal AGN 
contribution in the fitting regions. Typically, we 
obtain 68% confidence intervals with a width of 
10-20 km s-i. 
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Fig. 3a. — Best-fit models (red trace) overplotted on the host galaxy spectra (black trace) in the two spectral 
regions described in the text. Top of each panel: Rest-frame host spectrum with a fitted broadened stellar 
template overplotted. Bottom: residual from the fit (blue trace). 
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Fig. 3b. — Best-fit models (red trace) overplotted on tire host galaxy spectra (black trace) in the two spectral 
regions described in the text. Top of each panel: Rest-frame host spectrum with a fitted broadened stellar 
template overplotted. Bottom: residual from the fit (blue trace). 
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In reality, the uncertainty in the velocity disper- 
sion is larger than this because not aU features in 
the stellar templates used match the host galaxy 
features equally well. In order to estimate the un- 
certainties related to template mismatches we fit- 
ted a range of stellar templates to the host galax- 
ies and obtained a standard deviation of the mea- 
sured values of 0-* . We list the standard deviations 
for each range of fits for each region in Table [H 
As mentioned above, the systematic uncertainties 
due to template mismatches are generally much 
smaller in the red region than in the blue region, 
in agreement with the results of Hiner et al. (in 
preparation). 

The velocity dispersions measured from the 
blue and red regions agree within the uncertain- 
ties. In order to get representative velocity disper- 
sions for each object, we average the two values 
weighted by their uncertainties and the goodness 
of fit (as given by xled)- The weighted averages 
with their uncertainties (including systematic un- 
certainties added in quadrature) are listed in the 
last column of Table [H 

3.2. Galaxy Morphologies 

We analyzed the HST imaging data with 
the two- dimensional modeling program GALFIT 
(jPeng et a l. 2002) to achieve three primary goals: 
(1) to obtain photometry of the nucleus to be used 
in the measurement of A^bh, (2) to determine the 
size and characteristics of the bulge component, if 
present, and (3) to search for signs of interaction 
in order to properly interpret our results. 

It is imperative that the appropriate PSF is 
used when modeling galaxy morphologies. Em- 
pirical PSFs are ideal for this purpose, as they re- 
flect the same instrumental conditions that were 
used for the observations. We searched the HST 
archive for high signal-to-noise stars (typical peak 
counts were ^ 2100) that were observed with the 
same instrument and filter and that fell near the 
same CCD position as the QSOs. The input PSF 
image must be large enough to encompass the 
wings of the function, so we chose a large region 
and interpolated across extraneous sources within 
the region. Since the PSF is undersampled in the 
PC image, we broadened it with a Gaussian of 
a = 0.8 pixels. This increased the full-width at 
half-maximum (FWHM) of the PSF to be greater 
than 2 pixels, while preserving the encompassed 



flux. 

The science images need to be large enough 
to accurately estimate the background sky level. 
This is especially important, because the edges of 
galaxy profiles are sensitive to the background sky 
level, and can affect, for example, the Sersic index, 
n, of the profile fit to the galaxy. We cut each 
image to 15". 75 by 15". 75 (350 pixels x 0.045 arc- 
sec pix~^) centered on the science target, which is 
several times the size of the targets themselves. 
This large region sometimes included other ob- 
jects, which we either masked from the fit or fit 
with an extra component. To match the resolu- 
tion of the PSF, we also broadened the science 
images with the same Gaussian used on the PSF. 

We followed an iterative process while fitting 
the morphologies of the galaxies. Each fit included 
the point source (a scaled PSF) and background 
sky level estimated by GALFIT. We began with 
simple models of single bulge or disk components 
for the host galaxies and added complexity as nec- 
essary. All but one of the galaxies required two or 
more components to model the host. Where the 
image showed multiple objects in the field, the sec- 
ondary object was fit with its own component or 
masked from the fit. We generally used free index 
Sersic profiles to model the bulges, but in many 
cases, GALFIT did not converge on a reasonable 
Sersic index. To fit these hosts, we ran several it- 
erations of GALFIT holding the index constant at 
values n = 1, 2, 3, and 4. In each case the n — A 
fit the image the best. We adopted the best-fit 
models based on the values and residuals of 
the model subtracted image. 

We report magnitudes for the best fits of each 
of the targets. Visual inspection of the resid- 
ual images indicates that the point source and 
PSF were slightly mismatched. In order to esti- 
mate the uncertainty in the magnitudes, we also 
performed a direct PSF subtraction and deter- 
mined the upper and lower limits for which the 
PSF was clearly over- and under-subtracted, re- 
spectively. The magnitudes measured by direct 
subtraction were, in each case, within 0.1 mag 
of the value obtained from the GALFIT fits, ex- 
cept for 005055.7-^293328, for which the differ- 
ence was 0.36 mag. Our measured magnitudes 
are systematically f ainter than those reported by 



Marble et al.l ()2003() . with the differences ranging 



from less than 0.1 mag to 1.1 mag in the most 
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Table 3 
Best-fit velocity dispersions 



QSO 






Blue region 






Red rej 


;ion 


Weighted 


(2MASSi J) 


(T 


* 


Xred 


std. dev. 


a 


* 


^red 


std. dev. 


avg. cr* 


005055.7+293328 


169 


+6 
-4 


1.71 


15.4 


179 


+ 7 
-15 


0.69 


18.3 


176+}4 


015721.0+171248 


250 


+ 12 
-8 


6.75 


53.0 










250^^1 


022150.6+132741 


115 


+9 

-7 


1.96 


40.7 


140 


+7 

-7 


3.31 


20.0 


129^^2 


034857.6+125547 


136 


+9 
-17 


6.16 


26.1 


157 


+9 
-8 


2.99 


17.9 


152+1^ 


163736.5+254302 


159 


+2 
-8 


1.49 


16.7 


124 


+ 14 
-11 


0.51 


8.4 


12911^ 


165939.7+183436 


188 


+ 13 
-22 


0.50 


22.1 


150 


+32 
-26 


0.20 


16.8 




225902.5+124646 


198 


+ 10 
-2 


0.87 


25.8 


190 


+ 12 
-12 


0.69 


27.8 


194l^J 


230442.4+270616 


165 


+3 
-4 


1.53 


35.3 


125tl 


0.76 


3.1 


127+^ 


232745.6+162434 


229 


+7 
-6 


2.23 


21.6 


195 


+11 
-12 


0.89 


23.1 





Note. — Velocity dispersion in units of kms^ obtained for the blue and red fitting 
regions. The errors for the blue and red regions are 68% confidence intervals; Xred 
reduced for the best fit. The standard deviation in km s~^ of a series of measurements 
using different stellar templates is given for each region; this is an estimate of the system- 
atic uncertainty due; to template mismatches. The last column is the average of the cr* 
from the blue and red regions, weighted by their corresponding uncertainties. The errors 
given for the average values include the standard deviations added in quadrature. The 
blank spaces denote regions where a good fit was not obtainable. Note that the velocity 
dispersions in this table have not been aperture corrected. 
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extreme case (165939.7+183436). The differences 
likely arise from the fact that we fit simultane- 
ously a PSF with the different components of the 
host galaxy (e.g., bulge and disk), so the nuclear 
fluxes that we measure are less likely to be con- 
taminated by light from the host galaxy. Another 
reason for the differences in magnitudes may be 
the choice of PSFs: Whereas we selected empirica l 

(|2n03h 



Marble et al 



PSFs from the HST archive, 
used one of the QSOs in their sample as the PSF 
to subtract from all the other QSOs. They chose 
222202.2+195231 since it did not show significant 
extended emission after subtracting an artificial 
PSF created with the software package Tiny Tim 
(iKrist and Hookll200lh . 

In Table |4] we list the magnitudes of the PSF 
and Sersic components used to fit each object. 
The magnitudes are apparent HST F814W mag- 
nitudes, not corrected for Galactic extinction nor 
K-corrected, to allo w for direct compari son with 
resuhs pubhshed by iMarble et al.l (|2003[ ). Every 
object appears to have a clear bulge component. 
In Figs. Hal |4b] and |4c] we show the image, the 
best fit model, and the residuals after subtraction 
of each of the targets. One-dimensional surface 
brightness profiles of each object, including all the 
components used to model them, are shown in 
Fig. [51 

The object 005055.7+293328 is hosted by a spi- 
ral galaxy and has a prominent dust lane in the 
central region. In order to obtain the best estimate 
for the nuclear fiux we fit the host galaxy with mul- 
tiple components so as to minimize the residuals 
after subtraction. The best fit included the spiral 
arms, disk, bulge, and an additional narrow Sersic 
profile (not included in Table HI to account for 
the asymmetric light distribution due to the dust 
lane. While fitting all these components allows us 
to obtain a good fit for the PSF, it is likely that 
we are underestimating the fiux and the effective 
radius of the bulge, since we are splitting it into 
two components. 

For 163736.5+254302 and 165939.7+183436, 
we use additional components to fit the apparent 
companions. These components are not listed in 
Table SI but they are shown in Fig. [SJ 



3.3. Black hole mass estimates 

We estimated virial masses for the BHs in the 
sa mple by using the sc aling relation for Ha given 
bv lGreene et~aLl (|2010l) : 



Mbh = (9.7 + 0.5) X lO*' 



L 



0.519±0.07 



5100A 



1044 



erg s 



/FWHMhc 
V103 kms-i 



2.06±0.06 



Greene et al.l use t he vi rial coefficient mea- 
sured bv lOnken et all (|2004l ). i.e., < / > = 5.5, 



which is a factor of 1.8 higher than the value as- 
sumed for a spherically sym metric broad line re- 



gion (e.g., iKaspi et al.ll2000l ). This value is also 



consistent with that found by IWoo et all (l2010h 
from a sample of rever b eratio n mapped AGN. Al- 
though iGraham et al. find a lower virial 
factor (i.e., 2.8), we adopt < / > = 5.5 since 
this allows us to compare our results directly with 
those of previous studies (see [JH). 

In the scaling relation for virial masses, the size 
of the broad line region is a function of the FWHM 
of the broad component of Ha, FWHMhq, and the 
AGN continuum luminosity at rest frame 5100 A, 
L^im- Here we describe our procedure to measure 
each of these two quantities. 

3.3.1. Measuring FWHMna 

In order to measure the width of the broad Ha 
line, we first had to do a detailed decomposition 
of the narrow Ha component and the narrow N II 
AA 6548, 6583 lines. We performed all the line fits 
using the specfit task |Kriss 1994) in the IRAF 
STSDA S package, and follow ing t he method out- 
lined b y lGlikman et al. I (l2007h and[G reene and Hd 
(l2004l) . as described below. 

First, we created a narrow line model by fitting 
the S II AA6716, 6731 doublet. The line profiles 
often showed blue wings. Thus, we fit each line 
using two Gaussians, one of which was allowed to 
have skew (or asymmetry) , although we found that 
skew was not needed in all cases. The continuum 
was fit simultaneously with a power law function. 
To ensure the same model fit for both lines, and to 
avoid unrealistic degeneracies, we constrained the 
model for each line of the doublet to be the same. 
This was done by matching the widths of each 
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Table 4 

Results of Modeling the QSO Host Galaxies Using GALFIT 



QSO 


Component 


'TT-F814W 


^eff 


r-eff 


Sersic Index 


(2MASSi J) 




(mag)" 


(") 


(kpc) 




(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


nCtKCiKK '7_lOOQQOQ 

UUDUoo. / +zyoozo 


"DOT? 

sror 


on p;n 










bulge 


19.27 


0.22 


(J. 52 


4 (nxed) 




UlSK .'^ 




0.4D 


Q OO 


nit; 
U.iO 




arms 


io.OD 


1 CO 


4.oU 


U.Do 


Ui57zi.U+i7i24(5 


rbl" 


on A c 










bulge 


on 


u.iy 


U.D ( 


4 ^nxeuj 




disk 


iy.oo 


yj.i Z 


o 

2.0U 


i (nxeuj 


Uz2i5U.b+ioz74i 




21.17 










bulge 




2.oy 


1 .U4 


4 (^nxeuj 




disk 


1 Q CO 


O 1 Q 

z.io 


o.iy 


r\ t; 
U.o 


Uo4(30 ( .0+1/004/ 


T3CT? 

car 


01 C^T 

21.0 / 










bulge 


1 c\ on 

ly.su 


U. (y 


o 7n 


c; Q 

o.o 




companion 


20.88 


0.44 


1.51 


1 (fixed) 




companion 


20.08 


0.70 


2.40 


1 (fixed) 


163736.5+254302 


PSF 


21.37 










bulge 


19.88 


0.22 


0.92 


4.7 




disk 


19.87 


0.49 


2.04 


0.6 


165939.7+183436 


PSF 


19.11 










bulge 


18.21 


0.26 


0.76 


4 (fixed) 


225902.5+124646 


PSF 


19.40 










bulge 


18.82 


0.41 


1.33 


3.7 




disk 


19.51 


1.20 


3.89 


1 (fixed) 


230442.4+270616 


PSF 


20.70 










bulge 


19.23 


0.98 


3.65 


4 (fixed) 




disk 


21.50 


0.32 


1.18 


1 (fixed) 


232745.6+162434 


PSF 


22.12 










bulge 


19.20 


1.30 


6.56 


4 (fixed) 




disk 


20.07 




.',.22 


1 (Iixi-d) 



Note. — " F814W magnitudes not corrected for Galactic reddening or for 
intrinsic reddening. 
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Fig. 4a. — HST F814W images of the 2MASS QSOs. For eaeh row, the left panel is the original image, the 
central panel is the best model obtained with GALFIT, and the right panel is the residuals after subtracting 
the best model from the observed image. Each image is 10"xlO". 
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Fig. 4b. — Same as Fig. |4a| 
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232745.6+162434 



4c. — Same as Fig 



component and the relative positions of the dou- 
blet. Because we used a two-Gaussian model, 
an extra constraint on the flux ratio was neces- 
sary to ensure th e same total model for each line. 
iGarstanj (1952) calculates the transition proba- 
bility and intensity ratio of the doublet ''S - 
transition in S II. He finds it can range from 0.5 
to 0.8 depending on the parameters used in the 
calculation. However, he also cites a reasonable 
observed value of 0.74 based on NGC 7662. We 
adopted this observed value for our S II fits, as it 
produced the most reasonable fit of the range. 

In some cases the S II doublet was contami- 
nated by the Ha broad emission or telluric ab- 
sorption. In these cases we fit the [O III] AA4959, 
5007 doublet to create a narrow line model. Sim- 
ilar techniques as with the S II doublet were 
used to fit the [O III] doublet. In this case 
the doublet flux ratio was c onstrained to 2.87 
( Osterbrock and Ferlandl[20Q6 ). which was a good 
fit for our data. 

Wc then fit the region around Ha using the nar- 
row line model for N II A6548 and N II A6583. We 
constrained the relative position of the N II lines 
using the laboratory wavele ngths for these lines. 
While iGhkman et"al] (|2007t ) constrained the N II 
flux ratio to a theoretical value of 2.96, we found 
better fits when we let the flux of each N II line 
be freely fit. This often resulted in a flux ratio 
greater than 2.96. In addition, we included a nar- 
row line component of Ha using the fltted narrow 
line model and we used a power law to character- 
ize the continuum. Finally, we included a broad 
Gaussian (with initial width of 5000 km s~^) to 
model the broad line Ha component. We flrst 
flxed the broad Ha line to be symmetric, then 
added skew after its centroid position had been 
located by the fltting routine. We show the best 
flts with each of their components in Fig. [B] 

Specfit uses a simplex method to minimize . 
We adjusted input parameters and initial step 
sizes iteratively to probe different local minima 
and thus found the global minimum. The errors 
in the FWHM/f^ resulting from the flt are typi- 
cally ~15 km s~^. The true uncertainties, how- 
ever, are much larger since they are dominated by 
the uncertainties in the modeling of the narrow 
component and the decomposition of the broad 
and narrow components. 

In addition, although the AGN is much less ex- 



tinguished in these wavelength regions than in the 
regions where we measure cr*, there is still sig- 
nificant contribution from the host galaxy to the 
continuum around Ha. We found that subtracting 
different stellar Ha absorption lines (correspond- 
ing to different stellar templates or stellar pop- 
ulations) has no measurable effect on the value 
of FWHM^ffj. However, the shape of the stel- 
lar continuum around Ha is different for differ- 
ent stellar populations, and this can potentially 
affect the fitting of the broad component, espe- 
cially for the objects with the broadest emission 
lines. Thus the choice of stellar template to sub- 
tract from the AGN spectrum introduces a sys- 
tematic uncertainty. In order to estimate this un- 
certainty, we measured FWHMi/c from the differ- 
ent template-subtracted emission lines and found 
the standard deviation for the different measure- 
ments. We added this uncertainty in quadrature 
to the fitting and systematic errors. 

The final values for FWHMi^c, and their uncer- 
tainties are listed in Table [3 The FWHM//„ we 
measure for 225902.5-1-124646 is highly uncertain 
because Ha has a very asymmetric profile, with 
very extended emission on the blue side of the line. 
This is likely indicative of an outfiow, so that the 
A^BH we derive from it may not be reliable. We 
attempted a variety of ways to separate the core 
of the line from the outflow, by fitting the latter 
with multiple Gaussians blueshifted with respect 
to the core. In order to ensure that we obtain a 
conservative measure of jMbh, we chose the small- 
est value of FWHM/fa that would fit the line. The 
upper limit on the errors is given by the FWHM 
that would be obtained if a single Gaussian were 
used to fit the broad profile. 

We note that some studies, inclu ding those 



of re verberation mapped AGN (e.g., I Woo et al 



2010f ) use line dispersions rather than FWHM to 



calculate A^bh- However, these measurements are 
based on Hp rather than Ha. For our study, 
we are unable to measure the line dispersion or 
FWHM of H;3, since the AGN is still highly ob- 
scured at those wavelengths. Instead, we must 
use Ha as a proxy for H/3. The relations be- 
tween these two lines in the literature are currently 
much better constrained for FWHM than for line 
dispersions. For example, the rel ation between 
FWH Mff^ and FWHMj^^ used bv iGreene e\~^ 



(|201flf ) IS derived by fitting 162 objects 
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Fig. 5. — One-dimensional surface brightness profiles of the HST F814W images of the 2MASS QSOs, 
showing each of the components used to fit the objects with GALFIT. The black dashed lines in the residuals 
panels trace the one-sigma error bars in the observed data. 
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Fig. 6. — Continuum-subtracted, rest-frame Keck ESI spectra (black trace) of the 2MASS QSOs in the Ha 
region. The narrow line component models are plotted in blue, the broad emission lines in green, and the 
overall fit to the data in red. 
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fsee iGreene and Hoi l2005l for details), whereas 
McGill et all (|2008l ) derive their relation for Ha 
line dispersions from only 19 objects. Thus, we 
use the relation for FWHM//„ rather than that 
for line dispersions. This choice does not affect 
our main results: if we were to use line disper- 
sions rather than FWUMna, the resulting A^bh 
would be, on average, 0.08 dex greater. 

3.3.2. Measuring L5100 

The AGN continuum flux at rest frame 5100 A 
cannot be measured directly from our ESI spectra 
for two primary reasons: (1) The continua of red 
QSOs suffer from significant extinction at short 
wavelengths, and this extinction needs to be prop- 
erly measured. (2) The spectrum in this wave- 
length region has a large contribution from the 
host galaxy that also needs to be accounted for. 

In order to obtain the unobscured QSO flux 
at 5100 A, we first measured the extinction for 
each object by modeling its spectrum. As de- 
scribed in §3.11 we modeled the QSO continuum 
by including a polynomial (in x — ln{X) space) in 
the fitting procedure to obtain cr*. The resulting 
polynomial, C{x) x P{x), is then a good repre- 
sentation of the reddened QSO. In order to es- 
timate E(i3 — T^), we compared the polynomial 
to the Sloan Digitized Sky Survey ISDSS) com- 
posite QSO spectrum (jVanden Berk et al.l boOlh 
reddened with a Small Ma gellanic Cloud (SMC) 
redde ning law ( Prevot et al . 1984; iBouchet et al 



19851 ). with E(i3 — y) as a free parameter. The 
values oiFi{B — V) that we obtained using this pro- 
cedure are listed in Tabled Two thirds of our ob- 
jects span a range of E(B — V) values c omparable 
to those found bv lUrrutia et al. I (I2OO9I ) m a sam- 
ple select ed from the FIRST-2 MASS Red Quasar 
Surveys ijClikman et al.l 1200 ill , i.e., E(B - V) < 
1.5. The remaining three objects, however have 
E{B - y) - 2. Typical errors in E(B - V) derived 
from the fits arc '^O.l. 

Whenever possible, we also estimated F,{B — V) 
using Ha/H/3 ratios. Naturally, the values for 
these ratios become more uncertain with increas- 
ing extinction, so we only estimated F,{B — V) for 
the seven objects for which we could obtain reli- 
able ratios. The resulting F,{B — V) values were 
generally consistent with the values obtained from 
the fitting procedure, although the two could dif- 
fer by as much as E{B — y)^0.3. We fomid no 



trend in the way the two values differed; the av- 
erage F,{B — V) of the seven objects measured 
from the Balmer decrement was 1.49, compared 
to 1.45 as inferred from modeling the AGN con- 
tinuum. Since there is no a priori reason why the 
continuum and broad emission lines should suffer 
the same amount of extinction, we only used the 
Balmer decrement estimates as a rough compari- 
son to make sure the reddening estimates derived 
from the model were reasonable. 

As a further test, we subtracted the correspond- 
ing reddened SDSS composite spectrum from the 
spectrum of each of our targets to recover the host 
galaxy spectrum. There were some QSO features 
that were over- or under-subtracted in the result- 
ing host galaxy spectra. This is caused by the mis- 
matches between the SDSS composite spectrum 
and the intrinsic spectra of our QSOs. Indeed, the 
greatest source of uncertainty in the final values 
of isioo comes from such mismatches. However, 
the overall spectral energy distribution (SED) of 
each host galaxy matched the shape of the SED 
that would be expected for a galaxy with the given 
set of stellar absorption features observed in that 
galaxy. 

Once we had determined the intrinsic reddening 
for each object, we scaled each continuum model, 
C{x) X P{x), to match the observed fiux of the cor- 
responding QSO nucleus, as measured from HST 
WFPC2 F814W images using GALFIT (see Ta- 
ble Hj) and correct ed for Galactic extinction using 
the values given bv'Schle gel etaP (Il998l) . Thus we 
constructed continua of the reddened QSOs with- 
out galaxy contamination and taking into account 
flux loses due to the slit. Finally, we corrected 
these continua for intrinsic extinction (using the 
values listed in Table [5]) and measured the the 
flux of each QSO at rest-frame 5100 A. As men- 
tioned above, the greatest source of uncertainty 
comes from measuring F,(B ^ V) by making the 
assumption that the intrinsic continuum shape of 
our sources is equal to that of the SDSS composite 
QSO. The uncertainties for AL5100 listed in Table[5] 
include the uncertainty in the intrinsic power-law 
index, a\ (conserv a tively assumed to be '^0.05; 
Vanden Berk et aD ( 200l[ l) as well as the uncer- 
tainty in measuring E(i3 — V) by the process we 
outline above. 

The black hole masses derived from our mea- 
surements of FWHMifc. and AL5100 are listed in 
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Table 5 
Black Hole Masses 



Object 


QSO 


Weighted 


FWHM/fc 


E(B - V) 


AL5100 


Mm 


ID 


(2MASSi J) 


avg. (7* 


(km s"^) 




(10** erg s-i) 


(10** Mq) 


1 


005055.7+293328 


176+1*7 


2679ti75 


1.64 


2.52+1.24 


1 iq+O-40 


2 


015721.0+171248 


250tr4 


2196t^^ 


2.13 


18.74+2.05 


O/1+0.35 

Z.Zq:_g 72 


3 


022150.6+132741 


129ti 


4279tl^^ 


1.76 


1.78+0.24 


2 61+^ ''^ 

^•Di_i 01 


4 


034857.6+125547 


152+i^ 


4213tg5 


1.68 


4.16+0.95 


qo + 1.47 


5 


163736.5+254302 


129+1^ 


3178^^1^ 


1.93 


14.13+2.21 


4.15+1:?^ 


6 


165939.7+183436 




6742t?«f 


0.39 


0.94+0.24 


4.79t?:^3^ 


7 


225902.5+124646 






0.84 


4.03+0.55 




8 


230442.4+270616 


127+1 




0.99 


1.71+0.47 


r 77+1.62 
^- ' ' -2.02 


9 


232745.6+162434 


205+11 


2Ml+_f,\ 


1.93 


19.16+5.57 


2.59+_\ii 



Table m For completion, we also repeat the values 
of (7* from Table [3] in Table [S] As mentioned be- 
fore, using line dispersions rather than FWHM to 
estimate black hole masses results in masses that 
are, on average, 0.08 dex greater than those pre- 
sented in this table. 



4. Results 

Our results are summarized in Fig. [71 where 
the 2MASS red QSOs ar e plotted as black circles. 
We also plot objects from Woo et al. ( 2010[ ). which 
include the 24 local AGN that have published 
(7* and A^BH from reverberation mapping (green 
sq uares), as w ell as the A^bh-c* relation derived 
by I Woo et aL from that sample (green dashed 
line). For reference, the A^bh-' 7» relation for in- 



active local galaxies measured bv lMcConnell et al 



( 201l[ ) is plotted as a blue dotted line 

The majority of the red QSOs are clearly above 
the local relations. We have specifically com- 
pared our sample to AGN only to avoid any 
potential biases intro duced by compari ng active 
to quiescent galaxies (jLauer et al.l 120071 ) . and we 
have used a virial co e fficien t consistent with that 
found by IWoo et al.l (|20inh from the reverbera- 
tion mapped AGN sample. Woo et all (2010) es- 
timate an intrinsic scatter of 0.43 dex for their 
relation including all morphological types. In con- 
tras t, our objec ts are, on average, 0.80 dex above 
the Woo et aD relation, with six out of nine of 
them having log(A^BH) above the intrins ic scatter. 
Comparison with the iMcConnell et al.l relation 



for quiescent galaxies yields a similar, although 
smaller, discrepancy: our sample is, on average, 
0.65 dex above the relation, and with five objects 
being above the intrinsic scatter. 

When c ompa ring our objects to those of 



Woo et al 



(l2010h . it is pertinent to consider possi- 
ble differences in the width of Ha when measured 
from single epoch (SE) spectra versus rms spectra 



from reverberation mapped objects. iPark et al 
(|2012l ) find that H/3 is broader in rms spectra than 
in SE spectra, although the diff'erence between the 
two decreases as the intrinsic width of the broad 
lines increases 



Park et al] present a prescription 



to correct for this effect in relatively narrow-lined 
(<3000 km s~^ ) objects. To get an idea of how 
much this effect might impact our results, we can 
assume that the correction for H/3 holds for Ha, 
and that it is a reasonable approximation for the 
objects in our sample with broader Ha as well. 
W 'ith those as sumptions, if we apply equation (8) 
of lPark et all to correct the FWHlVIj^a in Table El 
we find that the black hole masses of the objects 
in our sample decrease by 0.13 dex on average. 
This correction, however, is not enough to bring 
the objects to agreement with the local A^bh-o"* 
relation, as th ey would st il l be, on average, 0.67 
dex above the IWoo et al.l (|2010t ) relation. Even 
if the correction for Ha and/or for objects with 
broader lines were greater, it is unlikely to be large 
enough to account for the offset we observe. 

It is important to consider other potential bi- 
ases that may be present in our sample. First, 
the morphology of the host galaxies could intro- 
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Fig. 7.— Red 2MASS QSOs (black circles) in 
the AAbh-o'* diagram. The green dash line is 
the Mbh-ct* for loc al AGN relation measured 
by IWoo et all (|2010l ). We also plot the Mbh- 
(Th. relation for quiescent l ocal galaxies measured 
by iMcConnell et all (|201l[ ) as a blue dotted line. 
Green squares are local AG N with reverbera - 

(I2OIOI) . 



tion mapped A^bh reported by I Woo et al 



Stars denote the Seyfert galaxies at z = 0.36 and 
z = 0.57 from the wor k by Woo et al. l2006t I2OO8 
and Treu et all ( 20071) . The two objects marked 
with hexagons do not follow the Faber- Jackson re- 
lation (see text for details) 



plying k-corrections for the spectral type of each 
galaxy and p assiye eyolution using th e on-lin e cal- 
culator from |yan Dokkum and Franx ( 2001 ) . The 
dashed line in the left pa nel of Fig. [51 is the 



Faber - Ja ckson relation from [Nigoche-Netro et aL. 
(|2010t ). iNigoche-Netro et alT " use a sample of 



~90,000 SDSS-DR7 early-type galaxies to calcu- 
late the Faber- Jackson parameters (i.e., A and 
B in log cr*= A - BMr) for different magnitude 
ranges. Our objects span a magnitude range of 
-20.6 < Mi < -23.4, or -19.8 < Mr < -22.6 
(AB magnitud es, using the trans formation equa- 
tions giyen by Jester et al. 20051). In th is mag- 
nitude range, iNigoche-Netro et al.l ()2010f ) find A 
= 0.152±0.00 9 and B = -0.976±0. 193. Note that 
the relation of Nigoche-Netro et al.l is measured in 
the SDSS r-filter, and we haye simply transformed 
it to Cousins-1, without accounting for potential 
wayelength-dependent differences in the slope. 



duce a bias. It has been shown that the precise 
form of the Mbh-<^* relation has a dependence 
on t he morphology of t he galaxies in question 
(e.g.. lGraham et al.luOllI) . and t hat pseudobulges 



do not follow the same relation (jKormendy et al 



20111 ). As can be seen in TablelH eyery host galaxy 
in our sample has a bulge component that is best 
fit by a Sersic of index n ~ 4. Although most of 
the galaxies also appear to haye a disk component, 
the majority (7/9) are clearly dominated by the 
bulge. In order to inyestigate how the presence of 
the disk may affect our measurement of ct, in the 
bulges, we teste d whether they follow the Faber- 
Jackson relation ( Faber and Jacksonlll976 ). In the 
left panel of Fig.|Sl we plot yelocity dispersions ys. 
absolute Cousins 1 magnitudes of the bulge com- 
ponents of the galaxies. The latter were obtained 
using the output magnitudes from GALFIT giyen 
in Table S] corrected for Galactic reddening, ap- 



Fig. 8. — Left: Faber-Jackson relation for the 
bulge component of the 2MASS red QSO hosts. 
The errors in Mj arc smaller than the symbols 
(0.01—0.04 mag). The numbers indicate the ob- 
ject ID number, as listed in the yarious tables. 
The dash ed line is the Faber-Jackson relation mea- 
sured by iNigoche-Netro et al. I (I2OIOI) from a sam- 
ple of ^^90,000 SDSS early-type galaxies in the 
luminosity range spanned by our sample in the 
SDSS r— band. Note that we haye not accounted 
for potential wayelength-dependent differences in 
the slope. Right: The bulges of the 2MASS red 
QSOs plotted on the Fundamental Plane. The 
solid line is th e orthogonal fit to the Fu ndamen- 
tal Plane from Hyde and Bernard! ( 20091 ) deriyed 
from -50,000 SDSS-DR6 early-type galaxies. The 
dashed and dotted lines correspond to the yalues 
that enclose 68% and 95% of the SDSS galaxies, 
respectiyely. 
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Two host galaxies fall clearly off the relation: 
those of 005055.7+293328 and 015721.0+171248 
(objects # 1 and 2, respectively), both of which 
have disks with luminosities comparable to, or 
greater than, their corresponding bulges. As de- 
scribed in §[33 the host galaxy of 005055.7+293328 
has prominent nearly edge-on disk with spiral 
arms and a dust lane that goes through the bulge. 
To mimic the dust lane, we fit this object with 
two separate components, so we almost certainly 
underestimat e the flux contained in the bulge. 
Marble et al ] (|2003h fit a single component to the 
host galaxy of this object and find a M/ = —22.3. 
This magnitude would bring the position of the 
object to agreement with the Faber- Jackson re- 
lation. The other object, 015721.0+171248, is 
the only one in the sample that has a disk that 
is more luminous than the bulge. The rotat- 
ing component from a dynamically cold disk can 
bias cr* toward higher values (when the disk is 
edge-on) or lower values (when t he disk is face- 

2011al) . Although the 



to Cousins-I, us i ng th e corrections prescribed by 
Rothberg etall (|2012h . The sohd line in Fig. E 



Bennert et al 



on; see e.g. 

disk of 015721.0+171248 has a lower inclination 
than that of 005055.7+293328, its velocity dis- 
persion may still be somewhat increased by the 
dynamically cold component. More importantly, 
015721.0+171248 is the only object for which there 
is significant QSO contamination in the stellar ab- 
sorption line spectrum, so that we were not able 
to obtain a reliable fit of cr* for the red region. We 
have marked both of these objects with hexagons 
in Figs. [7] and [TT] to indicate that the tr* mea- 
sured in these objects may not be representative 
of normal bulges. Interestingly, these are the two 
objects that fall closest to the local AGN A4bh- 
(T* relation, and removing them from the analysis 
only accentuates the offset of the red QSO sample 
with respect to the local relation. 

We also consider the position that the bulges in 
our sample occupy on the Fu ndamental Plane (FP; 
Djorgovski and David 1987 ). We obtain < ^ >e// 
as an output parameter from GALFIT, apply- 
ing k-corrections and passive evolution to 2 = 
as above, as well as correcting from cosmolog- 
ical dimming. We compare it to the F P rela- 
tion derived by iHvde and Bernardil (|2009l ) from a 
sample of ~50,000 SDSS-DR4 (with parameters 
from DR6) early-type galaxies at < z < 0.35. 
Here we use their orthogonal fit to the plane in 
the SDSS I— band and we have transformed it 



(right panel) corresponds to the coeffici ents given 
in Table 2 of iHvde and Bernardil (|2009l) . and the 
dashed and dotted lines correspond to the values 
that enclose 68% and 95% of their SDSS early- 
type galaxies, respectively. While, on average, our 
objects fall below the FP, they are consistent with 
the scatter in the relation. The two objects that 
are off the Faber-Jackson relation also appear to 
be the ones with the largest offset from the FP. 

Besides the effects of a potential dynamically 
cold component on cr,, we need to consider the 
possibility that we are measuring (7* in systems 
that are not dynamically relaxed. At least a third 
of the host galaxies in our sample show clear 
signs of interactions (Figs. Eil |4bl He]). Many 
of the numerical simulations that are success- 
ful at reproducing the TWbh-c* relation presup- 
pose tha t AGN activity is t riggered by merg- 



ers (e .g., iRobertson et al.ll2006a, b: Hopkins et al 



20061 ) ■ However, they predict the value (7* af- 



ter the system has reached dynamical equilibrium. 
To investigate the evolution of ct* during merg- 



Sticklev and Canalizol ( 2012f) use N-body sim 



ers, 

ulations of merging galaxies to perform a fiux- 
weighted measurement of cr, through a diffraction 
slit. They find that, typically, cr* tends to increase 
right after the first passage, followed by a period 
of oscillation about the value that cr* eventually 
reaches once the system comes to dynamical equi- 
librium. Thus, the observed value of cr, in our ob- 
jects could potentially be different from their cor- 
responding equilibrium value, depending on the 
precise merger stage at w hich we happen to be 
observing them. However, Sticklev and Canalizol 
estimate that, even in extreme cases, the observed 
value of cr, should fall between 70% and 200% of 
the quiescent value, and the probability that the 
observed value of cr* will actually be far from the 
equilibrium value is rather small. Therefore, even 
in the highly unlikely case that our entire sample 
is caught during an oscillation stage when cr* has 
a lower value, the predicted quiescent value would 
not be high enough to bring our sample to agree- 
ment with the local A^bh-c* relation. Moreover, 
using more sophisticated high-temporal resolution 
numerical simulations, Stickley and Canalizo (in 
preparation) find that the time at which the value 
of cr* shows significant deviations from the equi- 
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librium value preceeds the period of measurable 
accretion onto the black hole (i.e., the time when 
the object appears as a QSO). 

Another possibility is that dust in the host 
galaxy may be obscuring a fraction of the stars 
and this may result in a biased measurement of cr* . 
With the exception of 005055.7+293328, we find 
no clear evidence for significant amounts of dust in 
the host galaxies of the red QSOs. Their spectra 
show that, while the nucleus suffers from strong 
extinction, the host galaxies are not significantly 
reddened. As mentioned in § 13.3.21 after subtract- 
ing a reddened QSO spectrum, the overall SED of 
each host galaxy matches the shape of the SED 
that would be expected from the observed stellar 
absorption features. This indicates that there is no 
severe reddening in the host galaxies. Moreover, if 
the bulges of these galaxies were significantly red- 
dened, they would not fall on the Faber- Jackson or 
FP relations (see Fig.[5|). However, we cannot rule 
out the possibility that a modest amount of dust 
reddenning (E{B~V)<0.1) could be present in the 
hosts. Therefore, it is important to consider the 
effects that dust in the host galaxy, particularly if 
located close to the central regions of the bulge, 
would hav e on the me asurement of_av The sim- 
ulations of Sticklev and Canalizol ( 2012 ) include a 
toy model for dust consisting of a slab of gray at- 
tenuating material, and show that a lower value of 
(7* can be measured if the light from stars closer to 
the central regions is more extinguished than that 
of stars at large r radii. The specific toy model 
used by [Sticklcy and Canalizol leads to a decrease 
on the measured value of ct* <15%. Correcting for 
this factor would not bring the red 2MASS QSOs 
to the local TWeh-c* relation. While more realis- 
tic models for dust attenuation are needed to truly 
quantify the effect of dust on the measurement of 
(7*, we repeat that there is no evidence indicating 
that our host galaxies suffer from significant dust 
extinction. 

Finally, we consider the bias that may be in- 
troduced from the sample selection. The objects 
were selected based on their near-infrared colors. 
The selection criterion that objects must have 
J — Ks > 2.0 implies that their SEDs will be dom- 
inated by the AGN in the near-infrared; if they 
were dominated by the host galaxy, their colors 
would be bluer than this. Thus, the host galax- 
ies of the most reddened AGN could possibly be 



biased toward lower luminosity (and lower mass). 
We investigate this potential bias by looking at the 
offset of the objects as a function of near-infrared 
colors and reddening. First, we plot our objects 
on the black hole mass - bulge luminosity relation, 
A^BH- Lbuigc, in Fig.|9l where we have transformed 
the bulge magnitudes, M/ from Fig. [8] to as- 
suming an average R — I = 0.6 for the sample. 
The solid line in this fig ure is the loc a l rela tion 
for inactive galaxies from Bettoni et al. I (l2003l) . as 
reported by iKim et al.l (|2008l ) . Not surprisingly, 
our objects have overmassive BHs with respect to 
the local relation: they are, on average, 0.89 dex 
above the relation. Note that two of the objects 
that have the largest offsets with respect to the 
relation (objects 1 and 2) are also the two objects 
that do not follow the Faber- Jackson relation. If 
we do not include these two objects, then the sam- 
ple is, on average, 0.69 dex above the relation. 

If there were indeed a bias favoring less lumi- 
nous galaxies for more reddened objects, we should 
see a trend of an increasing offset from the 7V(bh- 
Lbuige relation, AMR^buige, with increasing J -Ks 
colors and with increasing reddening, E(i? — V). 
Figure [TU] shows no such trends and, in fact, some 
of the most reddened objects also have the small- 
est offsets in magnitude. Therefore, we do not see 
evidence for a bias introduced by the color selec- 
tion of our sample. 

5. Discussion 

Our sample appears to have a true offset with 
respect to the local A^bh-c* relation, in the sense 
that the black holes in red QSOs ar e overmas 
sive compar ed t o their host galax ies. I Woo et al 



(|2006il2008h and lTreu et all (|2007l) find similar off- 
sets for samples of Seyfert 1 galaxies at z ~ 0.36 
and z ~ 0.57; their results are plott ed along 
with ours in Fig. [71 IWoo et al.l investi- 
gate the A^BH-o'* relation at three different red- 
shift bins (<0.1, 0.36, and 0.57), and conclude 
that there is evolution at the 95% confidence 
level. This is at least in qualitative agreement 
with several studies of the A^bh- A^ bulge rela- 
tion at higher r edshifts. U s ing 51 qu asar hosts 
(mostly lensed), IPeng et ah (|2006al lbllcl) find evi- 
dence tha t the ratio A^b h /A^ bul ge increases above 
z = I. iBennert et all (|2011bl) find similar re- 
sults in a sample of 1 < z < 2 X-ray selected. 
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MR.bui,. (mag) 

Fig. 9.— Red 2MASS QSOs in the A^bh- Lbuigo 
diagram. The sohd hne is the local relation for 



inactive g alaxies fromlBettoni et al.l (|2003l ). as re- 
ported by iKim et aL ( 20081 ). The numbers indi- 
cate the object ID number, as listed in the various 
tables. The errors in Mr are, on average, 0.1 mag. 
The magnitude of object 1 is likely to be under- 
estimated due to the presence of a dust lane (see 
text for details). Note that objects 1 and 2 do not 
fall on the Faber- Jackson relation. 
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Fig. 10. — Offset from the A^bh- ^bui relation 
as a function of J — Kg color (left) and redden- 
ing (right). The numbers indicate the object ID 
number, as listed in the various tables. Note that 
objects 1 and 2 do not fall on the Faber- Jackson 
relation. 



broad-line AGN by deriving stellar masses of the 
bulge component from multi-filter surface photom- 
etry in HST images. The implication fro m these 
and other studies (e.g., Decarli et al. 2010f) is that 
at z ~ 2 host bulges were undermassive rela- 
tive to their BHs (compared to today) , suggesting 
that BHs grew first and bulges have been playing 
"catch-up." 

While our results appear to be in line with all 
these studies, they also present addit ional chal 
lenges to the evolutionary picture. I Woo et al 



( 20081 ) find a potential trend in the evolution of 
■A^BH-c*, if indeed the higher redshift samples 
are the direct progenitors of the local AGN. They 
show this clearly in their Fig. 3, where they plot 
the offset in .A/f rh with respect to the local, quies- 
cent sample of iTremaine et afl (2002). We repeat 
this plot, including our objects, in Fig. [TTl where 
we have calculated the offsets with respect to the 
local AGN A^BH -gH. relation from reverberation 
mapped objects ( Woo et al. 2010t )._ and we also 



include the l ocal AGN samp l es of iBennert et al 



Greene a nd Ho ('2006bl). The BH 



(|2011al ) and 

masses in the lGreene a nd Ho (2006b) sample have 
been increased by log(1.8) = 0.255 to account 
for the differen ce in the virial coefficient used by 
Greene and h3 compared to that used in the other 



studies. Thus, all the BH mass es in Fig.fTDassume 



a virial coefficient of / 5.5 (|Onken et al. 1 120041: 



Woo et a l. 2010). 



Taken at face value. Fig. [TT] would indicate that 
there is a steep increase in the offset between red- 
shifts of 0.1 and 0.2, and that most of the evolution 
has occurred in the last 3 Gyrs. Alternatively, in 
the context of models that fin d a non-causal origin 
for th e scaling relations fe.g.. Ijahnke and Macciol 
201 ih . most of the mass existing in the disks of 



these objects has been transferred to the bulge 
through mergers in the last 3 Gyrs. This would 
also be consistent with our result showing that the 
two objects with the most significant disk compo- 
nents are also the two objects with the largest off- 
sets in the AIbh- Lbuigc relation (Fig. [9]). However, 
the comparison between local and higher redshift 
AGN samples is not straightforward, as the latter 
have, on average, significantly higher A^bh than 
their local counterparts. Every object aX z > 0.2 
has A^BH > 10^ A^0, whereas there is a dearth of 
objects in this mass range at lower redshifts. In 
Fig. [m we plot aU the objects with A^bh > 10^ 
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Fig. 11. — Left Panel: A representation of the 
offset in Mbh with respect to the local AGN 
A^BH-o"* as a function of redshift, after IWoo at al 



( 20081) . The o ffsets are with respect to the 
Woo et al.l (|2010[) Mbh-ct* relation for local AGN. 



Open stars are the samples of Woo et al. ( 20061 
2008 ). crosses are the AGN from I Greene and Ho 
(|2006b), open squares a re reverberation mapped 
AGN (IWoo et al.M2010l). a nd ope n triangles are 
AGN from iBennert et al.l (|2011al) . The 2MASS 
QSOs (this work) are plotted as filled circles, 
bridging the gap between local samples and higher 
redshift samples. The two objects marked with 
hexagons do not follow the Faber- Jackson rela- 
tion in Fig. m Red symbols denote objects with 
log{MBii/MQ) > 8. Right Panel: Same as left 
panel, but plotting the offset in ct,. 



Mq with red symbols. If we consider the red ob- 
jects only, the trend with redshift becomes less 
significant. The red objects at z < 0.1 have an 
average offset of 0.40 dex, compared to 0.80 dex 
in our sample, and 0.66 dex and 0.61 dex, respec- 
tively, in t he Seyfert 1 samples at z = 0.36 and 
z = 0.57 of lWoo et al. l (l2006l . l2008h . 

The difference in A^bh range between the local 
and the higher- z samples is due, at least in part, 
to the fact that they are drawn from flux-limited 
samples. In these samples, less luminous objects 
are not detected at higher redshifts and, since 
AGN luminosity depends partially on Mbh, ob- 
jects with lower M bh are less likely to be picked. 
Lauer et al. (l2007h describe a potential bias that 



may be present in samples of AGN at the highest 
A^BH- They show that a bias may arise from the 
conspiratory nature of a steep luminosity function 
and intrinsic scatter in the AIbh-c* relation, lead- 
ing to a selection that favors higher mass BHs over 
lower mass BHs. Thus studies targeting luminous 
or massive AGN may find an offset from the lo- 
cal relation that is only due to this selection bias. 



On the other hand, iTreu et al.l (|2007n run Monte 
Carlo simulations to investigate the effects on 
a hypothetical selection function log( A^BH/Alf7)) 
> 7.9, whic h applies to the objects of IWoo et al 



(I2006ll2008h as well as our own. They compare the 
measured offset to a simulated input offset with re- 
spect to the local A^bh-o"* and find that the bias 
is negligible, unless the intrinsic scatter at higher 
redshifts increases dramatically (of the order of 1 
dex). However, the fact that we find an offset from 
the local relation in all of the objects with Mbb 
> 10* -Mq, both in the local and the non-local 
universe, may indicate that there is still some bias 
at play. 

In order to further assess potential biases, we 
examine the dependence of the A^bh-c* offset on 
AGN luminosities and Eddington ratios (Figs. [T^ 
and [T^ . For the reverberation mapped AGN of 



Woo et all (l2010l). we use th e values of Lsipp pub - 
lished bv IPark et al.l (|2012|). iBent z et al.' (l 2006l ). 
Kaspi et al.l (|2000t) . and lOennev et a l. (200^. For 
consistency, we transfo rm hna in the sample of 
Greene and H o' ("2006^ to Ls i no, us ing the rela- 
tions given by Greene and Hoi ( 20051 ) . The result- 



ing plot is shown in Fig. [TH We see a striking 
trend: The offset increases with luminosity, with 
virtually all the objects with log(L5ioo/erg s^^) 
> 43.6 being above the relation, including those 
with A^BH < 10* Af©- Clearly, the majority (two 
thirds) of high luminosity AGN are in the non- 
local universe, and this could be driving the trend. 
However, we note that all of the high-luminosity 
objects in the local universe follow this trend as 
well. 

One possible way to explain this trend might 
be if the radius-luminosity {Rblb. — L) relation 
somehow changed dramatically at high luminosi- 
ties, if, for example, a point is reached where 
the broad line region is matter-bound so that 
the radius-luminosity relation would flatten. A 
different relation might also be expected if the 
SEDs of high luminosity objects are different from 
those of low er luminosity object s. The data points 
in Fig. 5 of iBentz et al ] (|2009l ) suggest that, in- 
deed, some flattening of the Rblb. — L relation 
may occur at the highest luminosities. However, 
even if that were the case, the potential flatten- 
ing occurs at luminosities significantly higher than 
log(L5ioo/erg s~^) — 43.6. In fact, if we were 
to fit only the objects in the luminosity range 
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of 43.6 < log(L5ioo/erg s~^) < 45.3 (correspond- 
ing to the luminous objects of Fig. [121), '^^ would 
o btain a steeper rel ation, closer to the relation 
of iKaspi et al. I (l200d) . and resulting in somewhat 
higher black hole masses in this range. Thus, vari- 
ations in the Rblr — L relation cannot account for 
the trend we see in Fig. [T2l 

U sing semi-analyt ical models of galaxy evolu- 
tion, iPortinari et al.l 12011) find that more lumi- 
nous QSOs tend to trace over-massive BHs with 
respect to th e intrinsic relation (similar to the bias 
described bv lLauer"et al. 2007). However, for ob- 
jects at z < 0.5, their simulations predict that the 
parameters of QSO host galaxies closely follow the 
distribution of those of inactive galaxies, with very 
small offsets ('^O.l dex) in host galaxy mass. 

To measure Eddington ratios we obtain rough 
bolom etric luminosities by assuming Lboi = 10 
L5100 (|Woo and Urrvll2002[ ). The plot of ATWbh 
vs. Lboi/LEdd (Fig. [131) shows large scatter and no 
correlations. However, since Fig. [12] shows that 
there may be an offset in AAIbh at high AGN lu- 
minosities, it may be instructive to consider sepa- 
rately objects of lower luminosity (log (Lsioo/erg 
s^^) < 43.5; blue symbols in Fig. [T31) and those 
of higher luminosity (pink symbols). A trend ap- 
pears to surface: objects with low Lboi/LEdd have 
more positive offsets with respect to TVlBH-f* than 
those with higher Lboi/LEdd- We caution that the 
trends may be largely driven by the fact that we 
are dividing the objects in bins according AGN 
luminosity, e.g., a highly luminous AGN can have 
low A^BH only if it has a high accretion rate. 
Such trends of increasing Lboi/LEdd with decreas- 
ing A^BH have been obs erved in large flux-limited 
samples of AGN (e.g., iNetzer and Trakhtenbrot 
20071 ). However, we note that the trends we see 
in Fig. [13] are in the offset with respect to A4bh- 
(7, and not just in AIbh- 

A similar trend with respec t to TWBH-Lbuigc was 
observed bv lKim et al. ( 2008 ) in a study of 45 z < 
0.35 type 1 AGN. They find that, at a given bulge 
luminosity, objects with high Lboi/LEdd have lower 
■A^BH- They suggest that this may be explained 
in terms of the growth phase of BH: If Lboi/LEdd 
decreases as a function of time during a black 
hole growth episode, then BH with the highest 
Lboi/LEdd have more "catching up to do" than 
those with lower Lboi/LEdd- We searched for this 
trend by dividing the objects of Fig. [T3]in bins of 



stellar velocity dispersion (as a proxy for bulge lu- 
minosity; not shown in Fig. I13p . Given the large 
scatter, we found no significant trends or correla- 
tions. We also searched for this trend in our sam- 
ple by plotting AMfl_f,„;ge (from Fig. ^ against 
Lboi/LEdd- Once again, we found a large scatter 
and no significant trends. 

The combination of the potential biases in these 
samples with the trends that we have uncovered 
in this study prevents us from making any defini- 
tive statements on whether we see evolution with 
redshift in the A^bh-c* relation. Clearly, bet- 
ter statistics are needed to disentangle the depen- 
dence of A^bh-c* on mass, luminosity, accretion 
rates, and redshift. Studying samples with a wider 
range of these properties could help us to better 
understand the scatter and the different biases. 



For local AGN, Bennert et al.l (|2011al ) indicate a 
future study of ~75 additional AGN with A^bh 
> 10'' Mq that can be used to probe the local 
AIbh-c* relation at the highest A^bh and, likely, 
at high luminosities. In the non-local universe, we 
are conducting a study of additional red QSOs as 
well as other QSOs that are similarly suited for the 
study of A^bh-ct^, including objects in lower AIbh 
and lower luminosity ranges. We will present re- 
sults from o bjects at 0.2 < z < 1 in forthcoming 
papers fe.g.. lHiner et al.ll2012l ). 
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Fig. 12. — Offset from the A^bh-c* relation 
as a function of AGN luminosity at 5100 A. 
The 2MASS QSOs (this work) are plotted as 
filled c ircles. Open star s are t he non-local sam- 
ples of IWoo et ai] (l2006i 120081). cr osses are local 
AGN from lGreene and Hoi (|2006h[ ). o pen squares 



are lo cal reverberation mapped AGN (jWoo et al 
2010), and open tria ngles are local AGN from 
Bennert et al.l ( 2011a[ ). Red symbols denote ob- 
jects with log(A^BH/A^0) > 8. 



Fig. 13. — Offset from the A^bh-c* relation as a 
function of Eddington ratio. The 2MASS QSOs 
(this work) are plotted as filled c i rcles. Ope n 
stars are the samples of Woo et ahl ( 2006, 2008) , 
crosses are AGN from I Greene and Hg (2006b ). 
open squares ar e reverberation mapped AGN 
(IWoo et al.l 120101) ■ and op en triangles are AGN 



from lBennert et all ()2011a ). Blue symbols are ob 



jects with log (Lsioo/erg s ^) < 43.5 and pink 
symbols are the objects with higher luminosities. 
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